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ABSTRACT 

The properties of a partially ionized plasma in a long cylindrical tube subject to a 
uniform axial electric field are investigated. The plasma is maintained by an external 
ionizing source balanced by bulk and surface recombinations. Collisions between neutrals, 
whose density greatly exceeds the density of charged particles, and of neutrals with ions 
are sufficiently effective for their velocity distribution to be close to a Maxwellian with the 
same uniform temperature, independent of the external field. The behavior of the plasma 
is described by a collisional two-fiuid scheme with charge neutrality in the interior of the 
tube. Approximate nonlinear equations for the hydro dynamical moments are obtained 
from a Boltzmann equation in which electron-neutral, electron-ion and electron-electron 
collisions are all important. It is found that under certain circumstances the current, and 
the temperature of the electrons undergo a drastic change, with hysteresis, as the electric 
field is varied. 



PACS codes: 52.25.Fi; 52.65.Kj; 52.20.Fs; 05.20.Dd 
I. Introduction 

Instabilities are ubiquitous in strongly ionized plasmas. They dominate the behavior 
of such systems and their study forms the core of the subject. The origin of the instabilities 
lies in the nature of the plasma interactions: on the one hand they are long range and thus 
can produce strong cooperative effects and on the other hand they become 'weaker' locally 
at high energies (or temperatures) as manifested by the decrease of the Coulomb cross 
section with energy rise^ . The situation is different in weakly ionized cold plasmas, systems 
which have attracted much attention recently'^'"^. In such systems collective phenomena 
play a smaller role and instabilities are less common. Nevertheless there are cooperative 
phenomena in these systems too, which, as we have shown earlier'^'^ for homogeneous 
idealized systems and will show here for more realistic laboratory situations, can lead to 
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dramatic abrupt changes in the state of the plasma, when such systems are driven by 
external fields. 

There is a large literature on the behavior of plasma in a cylindrical geometry (see 
for example Refs.6,7,8,9). The reason for considering this system again here is that we are 
interested in a regime in which the electrons are colder than in the self-sustained discharge 
plasmas generally investigated. The plasma has a sufficiently high degree of ionization to 
make the electron- electron and electron-ion collisions important. The condition of low 
electron temperature permits us to neglect inelastic collisions but requires that we include 
volume recombination effects. 

The phenomenon in which we are most interested is a rapid, essentially abrupt, change 
in the electron temperature and current as the external electric field crosses a certain 
critical value which depends on the neutral-ion temperature, degree of ionization, etc. The 
phenomenon is related to the well known runaway effect in fully ionized plasmas caused 
by the decrease of the electron-ion collision cross section as the external field increases the 
electron energy-'^^'^-'^'^^. In the partially ionized gas the presence of the neutrals prevents 
such a runaway. In fact, if one neglects collisions between the electrons (e-e), as is done in 
the swarm approximation-^^, valid when the degree of ionization is sufficiently small, then 
the stationary distribution (in the absence of inelastic collisions) will be of the Druyvesteyn 
form^^ for which the current and electron mean energy are smooth functions of the external 
field. The situation is similar when the size of the plasma is smaller than the energy 
relaxation length of the electrons^. There are, however, other regimes, even in weakly 
ionized plasmas, where the effect of the e-e collisions is sufficient to keep the electron 
distribution close to a Maxwellian^'^'^. Such conditions can lead to cooperative abrupt 
changes in the temperature and drift velocity of electrons as the external field is varied. 
The origin of the phenomenon is made clearer if one starts with a kinetic approach rather 
than with a macroscopic description. 

While it is possible, even likely, that effects related to those discussed here have already 
been observed indirectly in the behavior of discharges it would be useful for both theoretical 

and practical reasons to have experiments in which the parameters can be controlled, so as 
to study the phenomena in a quantitative way. We expect that the transition will be seen 
as hysteresis in the current-voltage (I-V) characteristic of the plasma when the external 
electric field is slowly varying in time^'^. 

We shall consider an experimental arrangement consisting of a weakly ionized gas in 
a tube of radius R subjected to a constant external axial electric field E. The plasma 
is assumed for simplicity to be produced through uniform ionization inside the tube, by 
some external source, at a constant rate a. It is balanced by two kinds of recombination 
processes: a bulk one and a surface process on the tube wall. The main bulk recombination 
for the regime we are interested in, are^° three body processes involving two electrons 
plus an ion, and dissociative recombinations in which a metastable atom-ion complex 
recombines with an electron. The rates depend on the temperatures and densities of 
the electrons, ions and neutrals in a rather complicated way; see Ref.15. For the sake 
of simplicity we lump the two processes together and assume an effective rate of bulk 
recombination proportional to T~^/^, where T is the electron temperature. We ignore 
the dependence of this rate on the neutral and ion temperature and on the pressure 
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which we keep more or less constant. The recombination at the wall is also treated 
phenomenologically. In particular we assume that the energy is absorbed by the wall 
which is kept at a fixed temperature, see section 4. 

We have in mind here a situation in which the great majority of neutral atoms are 
some kind of noble gas to which may, or may not, be added a small amount of a more easily 
ionized second species, though we realize that in the latter case the analysis would be more 
complicated^^ . This will be reflected mainly in the rates of ionization a and recombination 
7 since we shall always consider a regime in which the density of electrons n(r, t) is much 
lower than that of the neutrals, N, but big enough, due to the great disparity between the 
electron mass, m, and the ion-neutral mass, M, for binary electron-electron collisions to 
dominate the energy exchanges in electron-ion and electron-neutral (e-i and e-n) collisions. 
This requires^ that 

am .,„s2 



^ AikTY « — « 1. (la) 

Here k is the Boltzmann constant, and a is the total electron-neutral particle collision 
cross section, which is taken to be a constant in our work. Putting in appropriate values 
for the parameters on the left side of (la) gives (see Ref.4), 

1 » n/N » 8 • 10"^, 4 • 10"^ 1.5 • 10"^ (16) 

for He, Ne, Ar plasmas respectively, when kT is approximately 1 eV and it decreases as 
for colder electrons. The upper bound relates to the fact that we ignore any collective 
self induced electrostatic or magnetic interactions. 

The ions in our model are assumed to have a uniform temperature, Tj, the same as the 
neutrals, while their density is ?i(r, t), i.e. the plasma is treated as locally quasi-neutral. We 
assume axial symmetry and longitudinal homogeneity so the spatial dependence is only in 
the radial variable r < R. The different mobilities of ions and electrons are compensated 
by an internally generated radial ambipolar electric field F(r, t). We are thinking of a 
quiescent, long positive column which fills the tube^. 

The reason for considering external rather than field induced ionization is that in our 
previous works^'^, in which wc considered a spatially homogeneous case with an a priori 
fixed plasma density, we found that as we varied the external field E, there was a transition 
in the electron distribution between regimes of weak and strong coupling to the ions. (Ref.5 
presents a rigorous proof of such a transition for a greatly simplified model system). This 
mechanism of a kinetic transition, which we investigate here in a more realistic physical 
situation, requires a low electron temperature and relatively weak electric field, which 
seems hard to achieve when the ionization is produced by the field, see Sec. 4. 

The basic idea in Ref.4 and here is to consider situations in which the collisions 
between electrons are strong enough to force their velocity distribution /(r, v,t) to stay 
close to a Maxwellian Mf with temperature T and drift velocity w. The values of T and 
w, which are simply related to the first two velocity moments of /, are then determined by 
self consistent "hydrodynamic" equations, i.e. we evaluate the integrals entering the time 
evolution of T and w with the help of replacing / by this Maxwellian. For the spatially 
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homogeneous case this yields ordinary differential equations in time for T and w which can 
be reduced to a couple of transcendental equations for stationary values T{E) and w(i?) 
yielding, in some cases, S-shaped curves as functions oi E = |E|. In the region where w 
decreases with E the system is unstable. The current and the electron temperature can 
therefore be expected to jump upwards (downwards) when the field increases (decreases). 
This will occur at different values of the field intensity, see Fig. 3: in the lower (upper) 
part of the loop the electrons will be supercooled (superheated) creating in this way a 
hysteresis loop. The origin of this behavior lies in a changeover from e-i to e-n coupling as 
the dominant factor when the electron energy is increased by the field. 

In the situation analyzed in this paper we obtain nonlinear partial differential 
equations for ?i(r, t), T(r, t), w(r, t) and F{r,t). Their stationary solutions are the main 
concern of this work. This is presented in section 3 following the mathematical formulation 
of the problem in section 2. A discussion of physical situations where the transition 
predicted by the two-fiuid model used in this paper might be observed is presented in 
section 4. This is a much more restrictive domain than that given by (1), it requires 
essentially n/N ~ 10"'' while N ~ 10^^cm~^ for light noble gases when Tj is close to 
room temperature. At the transition the electron temperature jumps within a range of 
several T^, staying significantly lower than the usual temperatures in self-sustaining gas 
discharges. This determines the experimental set up we study here. The nature of the 
various approximations made is discussed in section 4 and in the Appendix. 



II. Mathematical Description 

The behavior of the system is described by kinetic equations for the electron and ion 
distribution functions /(r, v, t) and (r, v, t) normalized to the same density n(r, t). Under 
the above assumptions the Boltzmann equation for electrons has the form"* 



dfjr, V, t) _ e_ 
dt m 



E + -F(r, t) ■ Vv/(r, v, t) + v Vr/(r, v, t) 



aip{v) - ^nf{r, v,t) + { {bn + /I) 



f{r,v,t) + 



kT^ df{r,v,t) 



mv 



dv 



hn ^ ?? — 

^Lfir, V, t) + y [/(r, V, t) - f{r, v, t)] + Q[f]. 



+ 



(2) 



v" I 

In (2) the first two terms on the right side represent the ionization and recombination 
respectively: V'(v) is the normalized distribution function of newly born electrons with 
mean kinetic energy mVo/2, 



V = Vn 



The third term with the prefactor m/M, represents the diffusion in the speed of the 
electrons due to collisions with ions and neutrals. The effectiveness of the e-i collisions is 
proportional to n and is strongly peaked at small speeds. The constant b is given'^'^^ by 
b = A'Ke^L/'m? (L ~ 10 is the Coulomb logarithm), while I = 1/A'"cr is the mean free path 
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of electrons in e-n collisions. The next two terms in the right hand side of (3) represent 
respectively the effects of e-i and e-n collisions on the angular parts of electron velocities 
with 

d d 
L = -Q^{c^6^,^ - c^Cy) — , (yu, = 1, 2, 3) 

where c = v — W and W is the ion drift velocity. We neglect the difference between |c| 
and |v|, in particular in the energy exchange term in (2), since |W| << v. For the e-n 
collisions we have defined 



as the sphericalizcd average of /(r, v, t). 

For simplicity wc have written the e-i and e-n collision terms when ions and neutrals 
have approximately the same mass M. In the case of using an easily ionized impurity these 
collision terms have to be modified in an obvious way. However, as long as the density of 
these impurities is of order of n, and n/N « 1, we do not expect them to substantially 
affect our result. Finally the last term Q[f] represents e-e collisions. We do not specify 
here the exact form of this term (although we have in mind a Landau form^*^'^^'^^), since 
it conserves the hydrodynamical moments corresponding to density, mean velocity and 
kinetic energy which we are going to derive from (3). 

Wc could but do not write down a similar kinetic equation for the ions. As already 
mentioned we assume them to have a distribution close to a Maxwellian with a uniform 
temperature Tj, equal to that of the neutrals. The assumption of charge neutrality further 
implies that both electrons and ions have the same density 



n{r,t) = J f{r,v,t)d^v = J Mr,vut)d% 



and the same radial drift velocity, as we shall see in the next section. The maintenance 
of this quasi-neutrality is the task of the radial ambipolar field F{r). The effect of 
recombination at the wall and other wall effects will be discussed in the next section. 



III. Hydrodynamic Description 

Many properties of the plasma can be expressed in terms of the first few velocity 
moments of / and /j. To obtain these we multiply (2) by 1, v, v'^ and integrate over v. This 
yields expressions for the time derivatives of the density n(r, t), drift velocity (v) = w(r, t) 
and {v'^) = 3kT{r,t)/m + w'^{r,t) which however involve unknown integrals over /. To 
obtain a closed set of hydrodynamical equations we follow the procedure in Refs.4,5, which 
is often used for coUisional plasmas, and compute all integrals over v with the help of the 
replacement 

f{r,v,t)^Mf{r,v,t) (4) 

where 

Mf{r,v,t) = n{r,t){m/27rkTf/^exp[-m{v - wf /2kT] (5) 
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is a local Maxwellian with the same r and t dependent first moments, n, w, (v^), as /(r, v, t). 
This assumption is motivated by the fact that the local Maxwellian Mf is the unique 
distribution which makes the e-e collision term Q[f] = and should thus be a good 
approximation (after some transient time) for situations in which Q[f] is the dominant 
term in determining the shape of the distribution function. This requires, at the minimum, 
that n/N satisfy (lb) and that E, a and 7 be small so that the energy exchanges in the 
e-e collisions are most important in determining the shape of the energy distribution of the 
electrons; the effectiveness of the e-n and e-i collisions being reduced by the factor m/M 
which is the physically small parameter in the problem. Eq. (5) can be rigorously justified, 
as was shown in Ref.5, for sufficiently idealized models. We believe that the evaluation of 
these integrals via Mf remains valid, at least approximately, in certain realistic situations, 
see section 4 and Appendix. 

Using (4) we obtain four equations for the electron density n, the longitudinal and 
radial components of w, Wz,Wr, the temperature T, and the radially directed electric 
field F. In writing down these equations we shall neglect higher order terms in the drift 
velocity w{r) as compared to the electron thermal velocity u{r) = \J?>kT /m. (This is a 
good approximation when m « M, see section 4 and Appendix). 

a) Electron density equation 



dt rdr 
b) Electron fiux in the z-direction 



dn d , . m? 

[rnwr) + a — c — g. 



(6) 



1 d{nw^) eE [6 bn + iu^/9l 

^7 \ = -2V -w^ 3 

n at m \ tt u-^ 

c) Electron flux in the radial direction 



cnWz 



(7) 



1 dinwj.) eF(r, t) 8/6 u cniVj. 

H = —- \l —Wr — 



n dt m 
d) The energy balance equation 

1 d{nu^ 



9 V TT ' Z 



15 2n 

— — (nw ). 



(8) 



+ n — [F{r,t)wr + Ew^] + — — {rnWrU^) 
ot m 6r or 



avf 



Here 



cn 
'2u 



6 m 

TT M 



^ \ bn + 8w^/9Z 

1- ^ —■ 



(9) 



u'^ (r, t) = 3kT{r, t) /m = (v^) - (r) , v'^ = 3kTi/m, 

and we have set 7 proportional to T~^/^, 7 = c/tt^, see Ref.6. The temperature 
independent parameter c describes qualitatively the rate of the volume recombination. 



6 



In the case of dissociative recombination c is typically "^^'"^^ in the range 4 x (10^^ — 
10^^ )cm^/sec^ and close to the lower value for pure Helium. The temperature dependence 
of 7 represents a combination of both dissociative process and three-body recombinations. 
The latter will be most important when the electron density n is large and T is low. It 
gives a reasonable rate in the vicinity of 300° K for the neutral and ion temperatures. 

A similar procedure with the ion kinetic equation yields (6) for the zeroth moment 
because charge neutrality implies the same n(r, t) and therefore the same Wr{r, t). The ion 
radial flux is described by the following equation 

1 d{nwr) eF{r,t) 8 j6 U cnwr 1 d . 2\ ^ ^ t 2\ 



n ot M 9 V TT <i u-^ 2>nor nr or 

Here U is the ion thermal velocity {MU^ = 3kTi) and li is the mean free path in i-n 
collisions with respect to momentum relaxation. Unlike in (8) we cannot, in (10), neglect w'^ 
compared with U^. We do however neglect the ionic contribution to the total longitudinal 
current and assume that the ions have the same temperature as the neutrals since they are 
in a much better thermal contact with the neutral particles than the electrons. We believe 
that in our range of electron temperatures we make only a small error in taking = const 
and neglecting the ion heating. Dependence of on r can be easily incorporated by adding 
one more differential equation similar to (9). 

The nonlinear partial differential Eqs.(6)-(10) have to be solved for the five unknown 
functions F{r),n{r),u{r),Wz{r),Wr{r). In the stationary state the set (6 - 10) can be 
reduced by the elimination of Wz{r) and F{r) with the help of (7), (8) and (10). We obtain 

2]/ 6 m gn{r) + 4u^{r)/9l ^ ' 



and three coupled equations 



■{rnwr) = a - c (12) 



rdr u 



m d f v?'+vf\ d 2\ (8 / 6m vi cn . 

^ o ] + Z:Jzi™r) = -\n\ -V7T + 7:R\^'>^r, (13) 



M dr \ 3 J rdr^ '"^ \9\TrMli 

u{r)n{r)-^ - i^)-^ = -^^ -nwr + A{r)/wr (14) 
for n{r),u{r),Wr{r). Here g = b+ \-\/^c and 

. 1 pK f eE\'^ u^ir)nir) a, o 5 o/ m 1 n^ir) 



2\l&\mJ ^n(r) + 4u4(r)/9/ 2'^ 3 3 u{r) 



6 m 

TT M 



,y2 1 

1 - 



V- , ,6n(r) + 8M^(r)/9Z 



■u^(r) 



^(0 • (15) 
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It is well known^'^'^ that in the hydrodynamic approximation the assumption of the 
charge neutrality leads to singularities near the plasma boundaries due to the space charge 
in the plasma sheath caused by the different thermal velocities of electron and ions and e-i 
recombinations near the wall. The singularity manifests itself in Eqs. (12-14) through the 
determinant of the linear system for the derivatives becoming zero when the 

radial drift velocity, at some r = r', reaches the value 



Wr 



5Te + 3T. 



3M 



(16) 



that is close to the adiabatic sound speed^. According to (8) the infinite derivatives 
imply an infinite electric field F{r) and hence a breakdown of the model. Following 
Refs.7,8 we interpret this to mean that charge neutrality cannot hold for r > r' and 
r' marks the boundary of the plasma sheath. On the other hand the thickness of the 
sheath is^ of the order of a few Debye lengths = ^JkTi/ Aiine^ ~ 10~^cm when 



Ti = 300K, n 



10 



12 



cm 



Assuming that R » Ad we can neglect the difference 



between r' and R in Eq.(16) and instead use (16) as the boundary condition for (12-14), 
as is proposed by Persson^ and other authors (see Refs.6,8). 

Our task now is to solve (12)-(14) subject to the conditions 



^(0) = 0, ^(0) = 0, Mo) = o 



(17) 



at r = and to Wr{R) given by (16) at r = i?. The functions n(r), u{r) are finite on the 
tube axis, r = 0, Eqs.(12), (14), (17) therefore imply 



dwr 



(0) = 



1 



[a-cn\0)/u^{0)] 



dr 2n(0) 

and ^(0) = 0. The last relationship is very important and we rewrite it as 



(18) 



_ gn{0)+4u^{0)/9l ( 
u^{0) n(0) I' 



6 m 
ttM 



1 - 



W2(0) 



n(0) 



6n(0) + 8u4(0)/9Z 



—a 



^0^- V(0) 



_ 2 n^{0) 
~ 3^ u{0) 



(19) 



Solving the set (12-14) with conditions (16-19) we determine the density and 
temperature profiles as well as the average temperature T{E) and the total longitudinal 
current I{E): 



T{E) 



2m 
3kR^ 



R 



u^{r)rdr, I{E) = 27re 



I rn{r)wz{r) 
Jo 



dr. 



(20) 



When we studied in Ref.4 the spatially homogeneous problem without ionization and 
recombination, we had c = 0, a = 0, n = const, u = const, and (19) served as the equation 
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for the determination of or the electron temperature, as a function of E. We have a 
similar situation here when l/R « 1. In this case one can disregard the recombination 
on the tube walls and we have again a spatially homogeneous problem n = n(0), u — u{0). 
The only difference is that we have now, in virtue of (12), 



n 



Defining 



SV 6 Vac V ™ / 4 ' 



Vi Vf 



(19) can be rewritten in the dimensionless form 







(jJ + X 

x^ 



5/2 



6^ 



{x^ -l){uj + 2x^/^) + {x^ -Q)/2 



(22) 



It is obvious from (22) that (p is not monotone in x and therefore T{E) is not single- 
valued if the parameter ut is big enough. In a Helium plasma where M/m ~ 7000, we take 
c = 4 • lO^^cm^sec"^, ^ = 10 and find that uj should be larger than about 113 in order to 
have a transition corresponding to an S-shaped T{E) curve; see Fig.l, where co = 250. For 
smaller c the critical value of u is lower. 

In terms of the ion temperature and the degree of ionization ut can be represented as 

4 

^ e n 

U) — TlL- 



This is equal approximately 2 • lO^n/N for Argon with the ambient temperature of the 
background gas. We see that the transition requires n/N ~ 10~^ or higher. Real values 
for the electric field and current density j can be found with the help of following relations: 



™ ■ (as) v^--' — " 

Here I is in cm, E and j are in volts /cm and mA/cm^ respectively. 

The T = const Approximation 

As a first step in solving the nonuniform case (12-14) let us assume that the electron 
temperature is constant, u{r) — u{0). It seems reasonable to study such a simplified 
problem both as a guide for the more general case which we shall consider later and as an 
approximation which often yields reliable results for the positive plasma column^. 

The energy balance equation (14) cannot hold now for all r < R, but at the point 
r = it reduces to (19), and our task is to solve n(12,13) with the conditions (16-18). For 
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fixed E and n = n(0) Eq.(19) determines u and therefore the electron temperature in the 
tube. We then integrate (12,13), find the profiles n{r), Wr{r), compute Wz{r) with the 
help of (11), and using (20) find the total current I{E), which gives the current-voltage 
characteristic of the plasma. We can neglect the difference between b and g, because in 
noble gases g/b — 1 ~ /x ~ 10~^. 

To solve Eq.(12,13) numerically we pick a value of u and find a suitable n(0) which 
allows to satisfy (16). Eq. (19) is used for the calculation of E for each choice of u and 
n{0). In this way we obtain the functions n(r), Wz{r) and substitute them into (20). 
Fig. 2 shows the electron temperature and the electric current thus computed versus E'^ in 
relative units when l/R — 0.2. Curves T[E) are single- valued for a; < 66 with the same 
0, c as before in a Helium plasma. 

const 

We solved numerically Eqs. (12)-(14) with boundary conditions (16-19). Technically 
we choose for each ti(0) a trial value of n(0), compute E with the help of Eq.(19), and 
solve the differential equations (12-14) for this triple. The resulting Wr{r) will generally 
not satisfy the boundary condition (16) at r = i? and we then iterate with a different n(0). 
This search can be easily optimized, the procedure usually converges very fast and yields 
the profiles of n{r) and u{r) for r < R. From these we can obtain F(r), Wz{r), Wr{r) for 
a given external field E. In addition we also find the total current and the mean electron 
temperature over the tube cross section (20). The parts of the T{E) and I{E) curves 
with negative derivatives are unstable and physically non-accessible"^'^, which also shows 
up in the computation being extremely unstable there when we try to find the solutions 
numerically. Thus for E = E' in Fig. 3 and any trial temperature on the tube axis inside 
the gap region (0.046-0.065 V), T(r) immediately jumps to a point corresponding to T' 
after a few steps of integration. This difficulty does not occur when we keep T = const. 

A qualitative description of the results is as follows: 

a) n(r) always decreases monotonically. For some range of parameters a kinetic transition 
takes place, n(0) is then not a single-valued function of the external field or the mean 
electron temperature. The decrease of the volume recombination rate c enhances the 
transition. 

b) The I-V characteristic for the total current 1(E) and the mean electron temperature 
show transition like behavior for a; > 60. This is close to that obtained for the model 
T = const. 

c) T(r) is a smooth monotonically decreasing curve when E is large; for very small E the 
electron temperature passes a maximum at some r < R. This maximum is caused by the 
electron heating in the ambipolar radial field F which is equal to zero on the tube axis. 
At stronger fields E this effect is disguised by the external field producing more heat near 
the axis where the electron density is higher. 

d) The radial drift of electrons and ions as well as the ambipolar electric field F[r) rise 
strongly near the tube walls, the latter even goes to infinity in this model but the wall 
potential with respect to the tube axis is finite as one can see from Eq.(8). 
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We present a few curves for a; = 90 in Figs. 3, 4 which iUustrate the transition, 
the temperature profiles across the tube cross section and the radial dependence of the 
transverse electric field F{r). As E is varied we expect jumps in T{E) and I{E) at the 
arrows in Fig. 3. 

IV. Discussion 

In this work we investigated properties of the stationary state of a partially ionized 
plasma in a cylindrical tube maintained by a source of ionization and subjected to an 
external axial electric field E. Our calculations were done in a coUisional two-fiuid scheme. 
The neutrals are treated as a given background with fixed density N and temperature 
while the ions were assumed to have the same uniform temperature Tj and a density n(r) 
equal to that of the electrons, i.e. charge neutrality in the interior of the tube. 

The assumption that the heating of the electrons by the external and ambipolar fields 
together with the injection of relatively warm electrons in the process of ionization does 
not heat the neutrals and ions very much seems reasonable when one takes into account the 
density and heat conductivity of the neutrals and the strong coupling between them and 
the ions. To see this in a semiquantitative way we note that the temperature distribution 
for the neutrals is governed by the equation 

where -P(r) describes the heating sources and A is the thermal conductivity which can be 
taken to be constant. By keeping the walls of tube at a given temperature Tq we impose 
the boundary condition T{R) = Tq which gives 

1 f^dx '■^ 



T(r) = To + - / — / yPiy)dy. (24) 







Let us consider three cases: 1) P{r) is a constant, 2) linearly falls to zero, and 3) approaches 
zero parabolically. Then letting 



2 

P = ^ J P{r)rdr 



be the average intensity of the heat source, Eq.(24) yields, on the tube axis 

Ti(0) -To + l^, UO) = To + T3(0) -To + l^ 

respectively. The condition 

P « 2XTq/R^ (25) 

then guarantees an approximately uniform temperature of the neutrals. In fact we take 
To = Tj and Eq.(25) can be rewritten in the form 
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if, for simplicity, we neglect the effects of the neutral gas non-ideality and only take into 
account the heating by the external field, P = Ej. For a typical example of the situation 
considered in this paper, which will be described later in this section, Ej is four orders of 
magnitude smaller than the right hand side of (25). 

The assumption of a constant ion temperature will hold only approximately. We 
expect that away from the tube axis the ion temperature will become somewhat larger 
than the temperature of neutrals but not significantly so, except maybe in close vicinity 
of the wall^'^'^. Since the transition we are studying here takes place when T is only a few 
times higher than that of the neutrals, we may assume a constant ion temperature Tj, at 
least for the majority of ions. 

The kinetic model on which the hydrodynamic equations (6-10) are based is the 
nonlinear Boltzmann equation (2) for the electron distribution function /(r, v, t). To obtain 
the fluid equations we have used the approximation of replacing / by Mf in evaluating the 
integrals for the time derivatives of the hydrodynamic variables. The justification given 
in section 3 is the dominant effect of e-e collisions. But, as is clear from the behavior of 
the e-e collision cross section, which decreases for large speeds as f f{v) will not be 
uniformly close to Mf{v) for all v. For large speeds, the behavior of / in the stationary 
state will be dominated, for \E\ not very small, by collisions with neutrals which lead to a 
Dryuvesteyn type tail of the form^^ / ~ e~^'" 1^ . Very slow electrons, on the other hand, 
are strongly influenced by the e-i interactions and their distribution is distorted near v = 
also. However, since we are only interested in the evaluation of averages of quantities which 
are not too peaked at very slow or very high speeds, the main contribution will come from 
typical values of v, w ~ m = ^J?>kT /m, where the bulk of electrons are. In this region of 
velocity space / will be close to Mf whenever the "relaxation time" Tee associated with 
e-e collisions is small compared with other relaxation times in the problem which would 
try to drive the electron distribution away from a local Maxwellian. We now evaluate the 
different relaxation times for typical electrons in the spirit of plasma kinetic theory^^'^^. 

The frequency of e-e collisions-^ is bn/u^ so Tee, the time of energy and momentum 
relaxation in such collisions, can be approximated as 

Tee ^ 7-- (26) 

on 

The relaxation time of the newly created electrons, whose mean speed vq is higher than u, 
can be obtained from (26) by replacing u with vq — On the other hand an electron 

with speed v loses, in a collision with a slow ion having velocity Wj, a fraction 2m /M of its 
energy. Such collisions take place with frequency, v^i = bn/u^ = Tee~^, so the e-i energy 
relaxation time is 

M . . 

Tei ~ Tee- (27) 

m 

By the same reasoning, the energy relaxation time with neutrals, where the collision 
frequency z^en is v jl^ is 

Ten ~ M/mVen = Ml/mu. (28) 

We thus always have Tee/Tei = m/M « 1 and the requirement that Tee /Ten << 1 is 
just the left side inequality in (la). This assures that e-e collisions will play a dominant role 
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in determining the energy distribution of the electrons, at least for the case of a uniform 
plasma with fixed electron density n, considered in Refs.3-5. In the present case we also 
have to consider the electron lifetimes relative to recombination in the bulk, tb, and on 
the wall Tw- A simple calculation yields 

rs = — = -Tee ~ lOVee. (29) 

cn c 

To estimate tw we compare the flux of neutralizing ions to the walls, Uin{R) ■ 2ttR, and 
the total number of electrons in the tube cross section nnR^: tw ~ Rn/2uin{R). Taking 
roughly n{R)/n ~ 1/20 we obtain, after some substitutions, 

Tw ~ lOR/ui 10{R/l){MT/mTi)^Ten- (30) 

Thus for R/l > 100, > and the arguments used for fixed n should apply. 

We note here however that in order to get M ~ Mf with both T and w determined by 
the hydrodynamic equations, one actually needs to assume also that r^e is small compared 
to the momentum relaxation times, i.e. that (ree^^ei) ^ind (reei'en) ^ire small. This was 
assumed for the simplified model equations^ but clearly does not hold in the plasma. A 
better approximation scheme to solving (2) would therefore be to take / ~ Mj?, a local 
Maxwellian with mean drift w = 0. The leading corrections to this Maxwellian would be 
of order ^/m/M, and would have to be computed from a linearization of (2) about M°. 
The result would then be the input for the self-consistent evaluation of T. This is what 
was essentially done heuristically and approximately in Ref.4. To make the scheme work 
rigorously the effective strength of E and hence w must also be of this order. This is indeed 
what we find here, for the range of E in Figs. 1-3, with w/u ^ ^JmjM . In fact the results 
obtained in Ref.4 are qualitatively similar to those obtained with the present approximation 
scheme for uniform systems, which facilitates computation; see also Appendix. 

Realizability of Transition 

We shall consider now, within the scheme of section 3 and the even simpler relaxation 
time estimates discussed above, the physical conditions necessary for observing the kinetic 
transition manifested in the ^'-shaped curves in Figs. 1-3. Since the origin of this behavior 
is a changeover from e-i to e-n collisions as the dominant mechanism for the dissipation of 
the energy which electrons gain from the external field, it is necessary that at the initiation 
of the transition, i.e. at the bottom part of the 5'-curve, the electron temperature, T^, 
should be low enough for Uei/^en to be large. The requirement that this ratio, which 
can be written as {A^inL /9N){e^ / a) / {kTi,Y , be large compared to one and that we still 
have n/N « 1 imposes a strong constraint on the temperature of the electrons T^. 
Assuming a Helium plasma with the neutrals and ions at room temperature, Tj = 300° ii', 
R = 1cm, 9 = {vq/vI) = 10, corresponding to an energy of 1/4 eV for newly created 
electrons, a = 5 x 10~^^ cm^ (see Ref.l7) the ionization fraction at the center of the tube, 
corresponding to the bottom in Fig. 3 is, n/N ~ 4 ■ 10~^. We took l/R = 0.2 which 
implies N ^ lO^^cm"^ corresponding to NkTi = 10~^ Torr and n(0) ~ 4 • lO^^cm"^. The 
electric field in Fig. 3 goes from zero to about 0.2 V/cm and the total current does not 
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exceed AmA. For a Neon plasma a is smaller by a factor of 9 and so is the requirement on 
the ionization fraction. 

Whether conditions like these can be achieved experimentally is at the moment an 
open question. What is clear however is that if we have to rely on ionization caused by 
the external field E then Tb will be about 1 eV and the electrons will never be coupled 
strongly to the ions for n/N « 1. To have a chance of seeing a transition corresponding 
to a crossover from e-i to e-n coupling we would have to consider an almost fully ionized 
plasma in which the remaining neutrals would prevent the runaway effect caused by the 
external field E when only e-e and e-i collisions are considered, (see Ref.ll). 

To gain a better understanding of the dynamics involved in this non-equilibrium 
kinetic transition, let us analyze, in simple physical terms, the stationary state of our 
system in the presence of the electric field E. We assume as before that in each e-i 
or e-n collision the fraction of energy lost is 2m /M and the direction of motion of the 
colliding electron is randomized. Under the action of the force — eE an electron between 
two successive collisions changes its velocity from v to v — eEr/m, where the mean time 
of the free flight, is 

fhn v\ ^ 

'=[v^^i) ■ 

When V corresponds to the thermal velocity u the drift velocity of electrons obtained from 
(31) is 

(This almost coincides with Eq.(13) for w^)- In the stationary state the average energy 
gain between two collisions 

2m 

should be equal to the average energy loss {m/M)m{u^ — vf), 

We see from (33) that u{E) is not a monotone function when q = bnl/vf is larger than 
about 12. This gives in a simplified form the origin of the kinetic transition. For this 
analysis to be reasonable the energies of the electrons must not be spread out too much, 
i.e. the bulk of electrons must have their speeds close to u. The e-e collisions provide 
the nonlinear cooperative coupling which brings about this condition. In the cylindrical 
plasma confined in a tube which was studied here the critical value of q is found very close 
to 10 for plasma parameters used in Figs. 1-3. 

An S-shaped T{E) or in other words a non-monotone behavior of the function E(u) 
is caused by the decreasing term bn/u^ in (33). The recombination of electrons near 
the tube walls makes the density n a decreasing function of u, due to the growth of the 
charged particle mobility with temperature and therefore enhances the observation of our 
transition. The opposite role is played by the volume recombination whose rate decreases 
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with the electron temperature {u here) and leads to a rise of density with u. We thus have 
a competition between these processes in the cylindrical geometry and may anticipate that 
the smaller is the bulk recombination the easier would be an experiment. 
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Appendix: Nature of Approximations 

We rewrite now Eq.(2) in a dimensionless form in order to better see the mathematical 
nature of the approximations made in our work. Setting e = y^m/M, t' = etvi/l, v' = 
v/wj, r' = r/R, ew' = w/vi, 5 = l/R,E' = {el/emvf)E, F' = {el/€mvf)F, X = 
lyjacjvl^ n{r') = nir^^fcfaJi^^ and dropping primes we get 

dj^^ _ ^ _ ^ _ ^^^^ ^ 



V ov 



dt 



(4jn/9 + v^)[f + 



]_dl 
3v dv 



+ Q^Lf + v{f-f)+Q[f], (Al) 



where uj is defined in (21) and we have 

J i/jd^v = 1, J viljd% = 0, J v'^i^d^v = e, 



fd v = n{r), J vfd V — en{r)w{r), J v fd v~n{r)[x (r) + e w (r)]. 

We restrict ourselves now to the homogeneous case {R — oo) with a fixed electron 
density, n, without ionization - recombination. Eq.(Al) then takes the form 



a/ _^ V /- — — 
dr dv 



Zv dv 

In the stationary state the equations for the moments have the form 



E + ^ y ^{^+'^ $(v)d'\ = 0, 



£;w + ^$(0)- 



V 3 



(A3) 
(A4) 



where q = bnl/vf as before and we have set $ = //n so that 

J v^d^v = ew, J v'^^d^v 
If one uses the substitution (7), which now has the form. 



2,22 

X + e w . 



$(v) ^ (3/27r)^/2a;"^exp 



2x^ 



(v — ew)^ 



{A5) 



and neglects e^{w/x)'^ compared to unity Eqs.(A3,A4) transform into the set 



E 



2w 6 



X"^ \ TT 



{q + 4x79) = 0, 



{AQ) 
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Ew = J- {q + 8x^/9) = 0. {A7) 

V TT X"^ 

In the stationary case (A6,A7) can be solved for w{E) and x{E) yielding non-unique 
solutions, when q > 7.6 (see Ref.4). The solution of (A6,A7) for x can be found from the 
relation ^ 

TT X° 

Eq.(A8) is close to Eq.(33), which was obtained in a rough approximation without referring 
to the kinetic equation. 

It is natural to expand /(v, t) for the spatially homogeneous case in a series in e 

/(v,t) = 5^/,(v,t)6^-. (^9) 
i=o 

Substituting (A9) into (A2) gives a set of coupled equations for fj where higher components 
can be expressed through lower ones. We do not use this method in the present work, 
instead we have solved Eq.(Al) using the method of moments with e = 1.4 x 10~^, d = 
0.2, A ~ 0.1, and © = 10. In order to have a closed set of differential equations we have 
taken 

/-(^)/ 

which comes from the shifted Maxwellian. We expect the effect of our additional 
approximation to be small. 
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FIGURE CAPTIONS. 



Fig.l 

The dependence of x (representing the electrons thermal speed) and the current 
density on (j) (the electric field squared) for the spatially homogeneous case {R = oo) 
when u) — 250. The dimensionless units are defined in Eq.(21) in the text. 

Fig.2 

The electron temperature T and total current / • i? (in mA ■ cm) vs the electric field 
E ■ R {va volts) when T is assumed constant in the tube cross section, Tj = 300° and 
u = 90. 

Fig.3 

Plots of the mean electron temperature and total current when T ^ const and u = 90. 
(The same units as in Fig.2). The hysteresis loop is indicated by arrows. At values of the 
field between the end points of the loop, like E', the computation leads to values on the 
lower or upper branches of the loop, determined by how close the starting point is to one 
of them. 

Fig.4 

The radial profiles of the mean electron speed and ambipolar electric field F[r) ■ R 
(in volts) when a; = 90 for two regimes: the solid lines correspond to a small field 
{E-Rf=i OMV, f/Ti 2.5), the dotted lines to a larger field {E-R^ 0.07V, f/Ti 8). 
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